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Chapitre 2 

Résolution des équations à une 
variable: f(x) = 

2.1 Position du problème et définitions 

Soit / une fonction réelle et Dj son domaine de définition (Df C H). Nôtre problème a la formu- 
lation suivante : 

trouver x € Dj tel que f(x) = 

Montrer l'existence d'une racine, déterminer le nombre des racines éventuelles de l'équation f(x) = 

et le calcul approché de ces racines posent des problèmes qui ne sont pas résolus par une méthode 

générale. 

Nous introduisons pour la résolution de ce problème plusieurs algorithmes dont on étudiera les 

caractéristiques (dichotomie, méthodes d'approximations successives, 

Newton-Raphson, . . .). 

Définition. Si s est une racine, de / alors /(s) = 0. Si / est p fois différent iable et telle que: 

f(s) = 0,/'(s) = 0, /<"-%) = et f"(s) î 
alors s est dite une racine de multiplicité p. 

Localisation d'une racine 

Il faut commencer par localiser la racine, c'est à dire trouver un encadrement de cette racine. 

a. Méthode graphique 

En étudiant le tableau de variation et les propriétés de la fonction (continuité, monotonie, . . . , 
graphe) essayer de localiser ainsi tous les intervalles qui contiennent une solution de f(x) = 

b. Par application du théorème des Valeurs Intermédiaires 

Si / est continue sur [a, b] C H, et si f(a)f(b) < alors il existe au moins un point s £.]a,b[ telle 

que f(s) = 0. 

Si de plus / est strictement monotonne sur [a, b], alors cette racine est unique. 

2.2 Méthode de dichotomie ou des bissectrices 

Soit / une fonction continue sur l'intervalle fermé [a,b] telle que f(a)f(b) < 0, il existe au moins 
une racine s dans l'intervalle [a,b] telle que /(s) = 0. Si de plus / est strictement monotone dans 
l'intervalle [a, b], alors s est unique. Notons que la mesure de [a, b] est b — a 

Dès qu'on a déterminé a = (Xq et b = 6q tels que f(a)f(b) < 0, la méthode de dichotomie consiste 
à considérer le point cq = ^- et à calculer /(cq) (voir figure). 
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1. Résolution des équations à une variable: f(x) = 



• si /(co)/(a) < alors s £E [a, cq], on pose ai = q.q = a et b\ = cq 

• si /(cq)/(6) < alors s G [co, o], on pose ai = Co et oi = fro = b 
mesure de [a, cq] est égale à ^^ et c'est aussi égale à la mesure de [co, b]. 

• Si s € [ao,co] on considère ci = ai t * = QQ + C ° et on étudie le signe de /(ao)/(ci) et celui de 
/(co)/(ci)- Si /(a )/(ci) < alors s e [a ,ci] sinon s e [ci,c ] 

• Si s G [cq,6o] on considère ci = ai ^ 1 = co ^ ° et on étudie le signe de /(fro)/( c i) e t celui de 
/( c <))/(ci). Si /(6o)/(ci) < alors .5 (E [feo.ci] sinon s e [ci,c ] 

• Dans tout les cas de figure on a: mesure [aosCi] = [ci,Co] = [bo; c i] — ~^r 

ar récurrence, on construit deux suites (a n ) et (b n ). La suite (c n ) est définie par c n = Q "^ " . On 
remarque qu'à chaque opération on diminue de moitié la longueur de l'intervalle contenant la racine 

b-t 
2" 

l'erreur est au moins -■ On a donc l'algorithme suivant: 



s. Au bout de n opérations la longueur de l'intervalle devient -t^t. Le coefficient de réduction de 



aa + b _ a + b 
2 ~ 2 



Algorithme de dichotomie. 

On définit deux suites (a n ) et (b n ) de la manière suivante: 
a = a, b = b c 
si /(io)/(co) < ai = a , 6i = c ci = 

si /( h o)/( c o) < ai = c , 6i = 6 ci = 
Si a n , b n sont connus alors a n+ i, b n+ i sont définis comme suit: 

on pose: c n = a "^ n alors; 

si /(c„) = on pose a„ + i = b n+1 = c„ +i 

sinon : 

si /(c„) > on pose a„ +i = a„, 6„+i = c„ 

si /(c„) < on pose a„ +i = c„ , b n+1 = b n 

Cela étant schématisé dans la figure suivante: 



a + c ai + bi 



2 2 

bo + Cq ai + oi 



(a,f(a)) 



a=ao 



-f 



f(a)f(b) < 



b=bb 



(b,f(b)) 



Figure 2.1: Principe de la méthode de dichotomie. 
D'après le graphe on remarque que: 



b — a 



2 

b — a 



u n "n ' 
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On en conclut donc que la suite (c n ) ainsi formée converge vers la racine s de f(x) = 0. c n est une 
approximation de s, l'erreur e n au bout de la nième opération est donnée par e n = \c n — s\. On a 
la majoration de l'erreur suivante: 

b — a 

K -s\< -^ 

Pour avoir une approximation de s par c n avec une précision e, il suffit qu'on ait: 

b — a 

\C n — S\ < — - < 6 



et donc 



\n('- 



hr2 



Exemple 

Soit à résoudre par la méthode de dichotomie f{x) = e~ x — x = 0. / est une fonction continue sur 
1R en particulier sur [0, 1], de plus /(0)/(l) < et / est strictement décroissante sur [0, 1] donc / 
admet une racine unique s g]0, 1[. D'après l'algorithme de dichotomie, on a le tableau suivant: 



n 


«n 


b„ 


C-n 


Cn+l _ c n 








1 


0.5 


0.5 


1 


0.5 


1 


0.75 


0.25 


2 


0.5 


0.75 


0.625 


0.125 


3 


0.5 


0.625 


0.5625 


0.0625 


4 


0.5625 


0.625 


0.59375 


0.03125 


5 


0.5625 


0.59375 


0.578125 


0.015025 


6 


0.5625 


0.578125 


0.5703125 


0.0078125 


7 


0.5625 


0.5703125 


0.5664062 


0.0039063 


8 


0.5703125 


0.5664062 


0.5683593 


0.0019531 


9 


0.5664062 


0.5683593 


0.5673828 


0.0009765 


10 


0.5664062 


0.5673828 


0.5668945 


0.0004883 


11 


0.5668945 


0.5673828 


0.5671386 


0.0002441 


12 


0.5671386 


0.5673828 


0.5672607 


0.0001221 



L'erreur au bout de la 12 eme opération est inférieure ; 
à 10~ 15 près est donnée par 0.5671432904097811. 



0.000122. La valeur approchée de s 



Remarques 

a.) L'erreur ei2 = JC12 

b.) c n est une approximation de s, 

de Terreur. Pour n = 11 on a: \c\2 ■ 



2T5TT 



: 0.000122 



2 n +i ■ 
■ en 



c„| donne aussi une information sur l'ordre de grandeur 
: 0.0001221. 



2.3 Méthodes des approximations successives. 

2.3.1 Point fixe. 

Soit g une fonction de IR dans IR et xo G IR, on engendre la suite: 



X2 



ff(»l) 



g(x„ 
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y=g(x) 




1. Résolution des équations à une variable: f(x) = 
y=g(x) 



s X 2 X, x„ 

Escalier convergent 



g(s) 



y=g00 






/ y=X 






"~^ 


^ 








[A 


Ë, 








^-^_^ 










X 


1 

Bsc 


l 3 

ar 


JOI 


x 2 J 

converger 


t 




x x i x 2 

Escalier divergent 




x l x x 2 

Escargot divergent 



Figure 2.2: Principe de la méthode itérative. 

graphiquement cela donne: 

Dans les cas convergents ( avec g continue ) on voit qu'à la limite x = g(x) 

Définition: on appelle point fixe de la fonction g de H dans H un réel x qui vérifie g(x) = x 



2.3.2 Conditions d'existence 

Pour qu'un point fixe existe il faut que : 

a.) l'ensemble de définition de g (D g ) et l'ensemble image de g (g(D g )) vérifient 

g(D g ) c D g 



b- 


y- gOO 


/ y= x 




b- 


y = g< 





/ y=* 


/ % 






10 Cg 






La courbe C g devra être 
dans la partie hachurée 












x 




/ 




x 




a b 






a b 



Figure 2.3: Condition nécessaire d'existence de point fixe g(D g ) C D g 

La condition "a." à elle seule est insuffisante car on peut avoir une fonction qui vérifie bien 
g(D g ) C D g et qui ne possède pas de point fixe comme le montre le graphe de la fonction suivante: 



g est discontinue en a 
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b- , — , — y 
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/ C o 



y=x 



R g C D g et il n'y a pas de 
point fixe. 



Figure 2.4: <? discontinue en a, g(D g ) C D 5 n'est pas suffisante 



Donc comme seconde condition pour qu'un point fixe existe il faut que: 
b. g soit une fonction continue. 

questions. 

1. Est-ce que g a un point fixe ? 

2. Ce point fixe est-il unique ? 

3. La suite (#7i) nÉ w -> converge- 1- elle pour tout xq G D g ? 

Des éléments de réponse (conditions d'existence, d'unicité et de convergence) sont formulés dans 
les théorèmes qui suivent. 

Théorème 1: Soit g une fonction réelle de variable réelle dont le domaine D g est un intervalle 
fermé borné [a, b] et dont l'image g(D g ) vérifie g(D g ) C [a.b] = D g . 

On suppose que g est continue sur [a, 6]. Alors g a un point fixe sur [a,b] c'est à dire qu'il existe 
un point s G [a, b] tel que s = g(s) 

Preuve: Soit / une fonction définie à partir de g par: f(x) = x — g(x). 

f(a) = a — g (a) < car g(D g ) C [a, b] ( g(x) > a pour tout x G [a, b]). 

f(b) = b- g(b) > car ( g(D g ) C [a, b], g(x) < b pour tout x e D g ). 

Or / est une fonction continue ( comme différence de deux fonctions continues) d'après le théorème 

des valeurs intermédiaires il existe s € [a, b] tel que f(s) = =$- g(s) = s. 

Remarques. 

• Si D g est non borné ca ne marche plus. Prenons par exemple la fonction g(x) = x + 1 
définie sur D g = H qui n'est pas borné et qui vérifie néamoins g(D g ) C D g mais pour tout 
x G D g , x ^/= x + 1. 

m Les conditions imposées ci dessus n'interdisent pas l'existance de plusieurs points fixes. Si on 
souhaiterait avoir une seule solution on doit supposer que g ne varie pas trop rapidement. 

Une façon d'imposer que l'intersection du graphe de g et de y = x soit non vide est de supposer 
\ g' \< L avec L < 1. On arrive ainsi à la notion d'application contractante mais la dérivabilité 
n'est pas nécessaire. Dans ce cas on n'a plus besoin de supposer que D g borné et ceci est contenu 

dans le théorème suivant: 
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Théorème 2: soit g une fonction d'une variable réelle dont le domaine de définintion D g est 
un intervalle fermé ( pas forcement borné). Donc D g = [a,b] ou [a, +oo[ ou ] — oc. a] ou H. On 
suppose: 

1. que l'image de D g par g vérifie: g(D g ) C D g 

2. et qu'il existe une constante L telle que < L < 1 

pour tous ï,!/é D g : \ g{x) — g (y) \< L \ x — y \ 
Alors: 

a. g a un point fixe UNIQUE s, 

b. La suite x n définie par xq éE D g : x n +i — g{ x n) converge vers s quel que soit xq éE D g 

c. pour n = 1,2, ... on a e n = \x n — s\ < ^^z I x i ~ ^o I- 
e n définie l'erreur à l'itération n. 

Preuve: 

Existence du point fixe: 

Si D g = [a, b]; déjà vu Théorème 1. 

Supposons que D g = [a, +oo[. Si g(a) = a, a est un point fixe, il n'y a plus rien à montrer, sinon 

g(a) > a (g(D g ) C D g ). 

g(x) < \g(x)\ = \g(a) + g(x) - g(a)\ < \g(a)\ + \g(x) - g{a}\ 

< \g(a)\ + L\x — a\ 

< \g(a)\ + L\x\ + L\a\ 

< c-\-L\x\ ; c = \g(a)\ + L\a\ 

donc 

g(x) < c + L\x\ 

Soit x > y^z ( conilïie L < l,x > 0) 
g( x ) ^ c + ^ < -^(1 — L) -\- Lx = x 

on a donc g(a) > a et #(x) < x, g admet donc un point fixe entre a et x. 
Il faut prouver que x > a 

En effet: î, > kM±£M = IsMI+U^H + H = _| a | + \eM+M > _| a | + ^ = a 

Si D g =] — 00, a] on prend g(y) = g(—y) : g est continue (comme composée de fonctions con- 
tinues ). On est donc ramené au cas précédent. 

!(z) g(o) = a 
(ii) g(a) > a 
(iii)g(a) < a 

(i) g(a) = a plus rien à prouver, le point fixe existe. 

(ii) g(a) > a: on a toujours \g(x)\ < c + L\x\ ou c = |<?(a)| + L\a\ Soit x > -^z on a c ^ onc 
g(x) < c + L\x\ < x g(a) > a et g(x) < x x > a; il existe donc un point fixe entre a et x. 
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1. 3 Méthodes des approximations successives. 17 

(iii) g(a) < a. on a 

-\g(%)\ < \g(x)\ 

-\g(a) +g(x) - g(a)\ < \g(x)\ 

-\9{a)\ - \g(x) - g{a}\ < \g(x)\ 

-L\x-a\ - \g(a)\ < \g(x)\ 

—L\x\ — L\a\ — \g(a)\ < \g(x)\ =$- —(L\x\ + \c\) < g(x) 

Soit x < ^Er =$- x < a et — (L\x\ + |c|) < g{x)g(x) > —L\x\ — \c\ > x d'où l'existence d'un 
point fixe entre a et x. 

Dans les 3 cas on a montré l'existence d'un point fixe, reste à prouver son unicité : 
Unicité du point fixe: 

En effet: soient s\ et 82 deux points fixes, \g(si) — g(s2)\ < L\s\ — $2] ==$■ |si — S2I < L\s\ — sa| 

Si si 7^ 82 j|si — S2I 7^ =r- L > 1 ce qui est absurde. 

Montrons b.) Pour montrer la convergence de (x n ) e jw vers s, considérons : 

\x n - s\ = \g(x n -i) - g(s)\ < L\x n _i - s\ < L 2 \x n _ 2 - s\ < . . . < L n \x - s\ 

< L < 1 ; lirn T1 _ i ._|_ 00 L n = 0; ceci montre que lim n ^ +co x n = s 

c.) donne un encadrement de l'erreur \x n — s\ après n itérations. Pour prouver (c) on doit montrer 
que 

|a>*+l - Xk\ < Lk \ x i ~ x o\ 
En effet 

\x k+1 - x k \ = \f(xh) - /(aîA-i)l < L \ x k - Xk-i\ < L 2 \xk-i - x k _ 2 \ < L k \xi - Xo\ 
maintenant soient m,n tels que m > n; considérons 

\ x m x n\ \ x m x m—l "T S'm— 1 x m—2 < ^m— 2---- *^n,+l ~T~ ^n+1 *"7l| 

_ j^m x m— 1\ ' \ x m—l Xjn—2\ i ... ~r |2-n+l x n\ 

< i m "Vi - x \ + £ m_2 |ici - x \ + ... + L"|ki - as | 

< (L m - 1 + L m - 2 + ... + L")| 2;i -x | 

< L n (l + L + L 2 + ... + L m - n - v )\x l -x \ 

1 Tm—n 

< £ ^ ki -x | 

< Xi — Xn\ 

- i-L 1 0| 

Si on fait tendre m — > oo tout en gardant n fixe, on a l'encadrement de c). 

Remarque: 

1. pour n fixe, l'erreur | x n — s | est d'autant plus petite que L est proche de 0. 

2. pour n fixe, l'erreur diminue très lentement si L £» 1. 

Proposition: Soit y une fonction dérivable dans [a, b], si g' vérifie : max xe r w |</(x)| = £ < 1 alors 
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1. Résolution des équations à une variable: f{x) = 



g est contractante sur l'intervalle [a.b]. 

Preuve: d'après la formule des accroissements finis: g(x) — g(y) = g'(^)(x — y), £ E]x.y[ donc 
\g(x) — g{y)\ = \g' (£,)\\x — y\ < L\x — y\ avec < L < 1. 

Remarque: Il faut que max|</(:r)| < 1 et non j</(#)| < 1 comme le montre le contre exemple 

suivant: 

si g{x) = x -\- - pour x G [l,+oo[, g'{x) = 1 — \. On a pour tout x G [1. + oo[ \g'(x)\ < 1 mais 

m&x.\g' (x)\ = 1 pour x > 1 on a aussi \g(x) — g{y)\ = \x — y\g'{a) < \x — y\ avec a G [1, +oo[ mais 

g n'a pas de point fixe. 



Théorème 3: g une fonction définie sur [a, b] = D g on suppose: 

• g([a,b]) c [a, 6] 

• m.ax.\g f (x)\ = L < 1 pour tout x G [a, b]. 



alors on a la majoration de l'erreur e n 



d'itérations qui assure e n < e. est n > N avec: TV = oge lo £ ge " 



s\ suivante: e n < L n eQ. Le noiubre minimal 



Preuve: 



\g(x n -i) - g(s)\ < L\x n _j - s\ < L n \x -s\ = L n e 



Proposition: Soit s une solution de g(x) = x avec g' continue : 

• si |</(s)| < 1 alors il existe un intervalle [a, b) contenant s pour lequel la suite définie par 
xq ÉE [a, 6] et x n +i = g(x n ) converge vers s. 

• si |</(s)| > 1, pour tout #0 ^ s la suite définie par xq et x^^i = g(x n ) ne converge pas vers s. 

Exemple: 

Soit à résoudre e~ x — x = <=> e~ x = x. 

pour x G [0,1] = D g - g(x) = e~ x : g([0, 1]) = [e~\l] C [0.1]. g'(x) = -e~ x , max\g'{x}\ = 1 sur 

l'intervalle [0, 1] n'est pas inférieur à 1, donc il faut changer d'intervalle. 

Prenons par exemple D g = [h,log2] on a g{\^,log2\) = [2, e ~ 2 } G [i,/o<72], max|</(:c)] sur \-^ 1 log2] 

qui vaut L = e~2 qui est bien inférieur à 1. 

Partant de Xq = 0.5 on a le tableau des itérations suivantes: 



n 


*Tl 


g(x„) 


■'')( + ! x n\ 


x' 





0.50000 


0.606531 




0.567624 


1 


0.606531 


0.545239 


0.106531 


0.567299 


2 


0.545239 


0.579703 


0.061292 


0.567193 


3 


0.579703 


0.560065 


0.034464 


0.567159 


4 


0.560065 


0.571172 


0.019638 


0.567148 


5 


0.571172 


0.564863 


0.011107 


0.567145 


G 


0.564863 


0.568438 


0.006309 


0.567144 


7 


0.568738 


0.566410 


0.003875 


0.567144 


8 


0.566410 


0.567560 


0.002328 




9 


0.567560 


0.566907 


0.00115 




10 


0.566907 


0.567278 


0.000653 




11 


0.567278 


0.567067 


0.000371 




12 


0.567067 


0.567187 


0.000211 




13 


0.567187 


0.567199 


0.00012 




14 


0.507119 
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l'erreur lors de la 14 eme itération est inférieure à: jzz\xi - scq| = „ 3 e 93 4 69 0-106531 = 0.000247 



2.4 La méthode de Lagrange ( ou de la sécante ) (1736—1813) 

Le principe de la méthode est le suivant: on localise deux points A = (a, f(a)) et B = (b, f(b)) tels 
que f(a)f(b) < (voir figure), on trace le segment de droite (AB) joignant A et B. L'équation de 
cette droite est donnée par: 

V = ./ ( a + (x-a) 

b — a 

soit Ai{xi,$) le point d'intersection de la droite AB avec l'axe des abscisses, Ai vérifie l'équation 

de la droite AB, donc 




B(b,f(b)) 



Figure 2.5: Principe de la méthode de Lagrange. 

Considérons le point i?i(xi,/(;ci) et traçons le morceau de droite (BiB) joignant Bi et B. 
L'équation de la droite (BiB) est donnée par: 

x/W-/(«l) 



y = f( x i) + (x-xi)- 



b — x\ 



soit A2(%2, 0) le point d'intersection de la droite (B\B) avec l'axe des abscisses, A2 vérifie l'équation 
de la droite (BiB), donc 



352=351- f(xi)- 



b — x\ 



'f{b)-f(x 1 ) 
même opération répétée plusieurs fois permet d'obtenir une suite (x n ) telle que 

b-X n _ f(x„) 



Xn+\ = g{x„) = X„ - f(x„) 



f(b) ~ f(Xn) 



b—x tl 



La suite (x n ) ainsi obtenue est convergente. Dans le cas où la dérivée seconde f" a un signe 
constant: la suite (x n ) est, soit croissante majorée, soit décroissante minorée. 
En effet, deux cas se présentent: 
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1. Résolution des équations à une variable: f{x) = 





Blb.f'lbn 



B(b,f(b)) A(a,f(a)) 



Figure 2.6: Méthode de Lagrange avec j"{x) > et f(a) > figure (a), f(a) < figure (b). 

a.) /" > et f(a) > 0: Dans ce cas l'extrémité a est fixe et les approximations successives 
sont: 

x = b 

#rH-l = X n — f{x n ) f( Xi j_f(a) 

d'après la figure on remarque que: 

a < s < . . . < x n+ i < x n < . . . < x-i < Xq 

Dans ce cas la suite (x n ) est décroissante minorée par a et est donc convergente. 

b.) f" > et /(a) < 0: Dans ce cas l'extrémité b est fixe et les approximations successives 
sont: 

Xq = a 

x n+l — x n ~ J \ x n) f /m_ f'(x n ) 

d'après la figure on remarque que: 

Xq < xi < . . . < x n < x n+ i < . . . < s < b 

Dans ce cas de figure la suite (x n ) est croissante majorée par b et est donc convergente. 
Remarque: Si /" < 0, on se ramène au cas précédent en considérant —f(x) = 0. 



2.5 Accélération de la convergence 

2.5.1 S 2 d'Aitken 

Quelques fois la suite (x n ) converge vers s très lentement mais la suite Xn +*~ s converge vers une 
constante (avec \9\ < 1). Alors on peut considérer la suite (x f n ) définie par: 



^n+1 x n) 



pour tout n > 



x n+2 ~ 2z n+ i +x n 
(x' n ) converge plus rapidement que (x n ) , c'est ce qui est exprimé dans le théorème suivant: 

Théorème 4: Si x n -> s (avec x n ^ s, \fn £ IN) et si x '^ l _~ s -^(9^0 (\9\ < 1) alors la 
suite (x' n ) converge vers s et on a: lim^^+co ^ ' - ~ s = 
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Preuve: Soit e n = x n — s alors 

(x n+1 -s- (x n - s)) 2 



x„ s — x n s 



(x n+2 - s) - 2(x n+1 - s) + (x n - s) 

( € n+l — En) 



f -n+2 — 2e ra+ i + e n 

X n -\-\— S Cïi+1 



ou: 



soit 9 n = n+ ]_ s = ^±1. d'après cette définition on a donc Cn+2 — ^n+i e n+i — &n+l^n e m d'où 

*-.-«.- rM^S ,,^ = i. (e "" 1)2 



En(ôji+iS n — 26„ + 1) e„ (ô„ + iô„ — 26„ + 1) 



or lim„^ +00 (s^f^J^!, = 1 donc 



lim ?"__! = 

7Î— > + OQ e n 



La méthode du rf 2 d'Aitken est appliquée pour résoudre e x — x = 0, les résultats sont donnés 
dans le tableau précédent. On remarque que l'on a atteint une précision de 10~ 4 au bout de la 
2 eme itération, alors que la même précision n'est atteinte qu'au bout de la 13 eme itération par la 

méthode itérative. 

2.5.2 Ordre d'une méthode 

On a donné précédement une majoration de l'erreur e n < jzxl* c i — x o\- On donne maintenant une 

évaluation asymptotique de Terreur: 



Théorème 5: Soit g une application de [a, b] dans [a, b] on suppose que : 

• g est contractante, < L < 1 

• g' est continue sur [a,b] ( g de classe C 1 ) 

• g' de signe constant sur [a, b] 

alors si eo 7^ on a: Vrc ÉE IN* e n 7^ et lim n ^, +00 e " + = |</(s)| 

Preuve: 

e ^ <£> x ^ s 

O g(xo) ^= g{s) = s (g strictement montotone) 

O xi 7^ s 

<^> ei t^O 

<£> : : 

& e n t^O 

d'après la formule des accroissements finis on a: 

e n+1 = x n+1 - s = g(x n ) - g(s) = g (£, n )(x n - s) = g'{^ n )e n 

avec Ç n Çz]x n ,s[, donc à la limite on a: 

km = \g(s)\ 

n— *+oo e 
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Théorème 6: Soit g € C 1 ([a, 6]), on suppose qu'il existe s G [a, b] tel que g(s) = s et g'(s) = 0. 
Alors il existe un e > tel que Vieq €]s — 6, s + e[ l'itération x n +i = g(x n ) commençant par xq 
converge vers s . 
De plus, si g est C (]s — e, s + e[) et si e n ^ pour tout n G IN on a 

lim ^±i = V(,) (2.1) 

n— ►+oo e- z 

Preuve: Soit < L < 1, la fonction 5' est continue et est nulle en s, il existe un intervalle 
/ =]s — e, s + e[ telle que |</(#)| < £ pour tout x € I donc p est contractante. De plus \g(x) — s\ < 
L\x — s\ < Le < e ceci implique que /(/) C I et donc lim n ^ + . 30 x n = s. 

e n+ i = g(x n ) - g(s) = (x n - s)g'(s) + -(x n - »)V(£n) 

= -(ar„ - s) 2 g"(Ç n ) avec Ç n e]# n ,s[ 

de la dernière égalité on tire: ^^ = %ff' f {£n) e "t à la limite lorsque n tend vers +00: Ç n tend vers 



s, on a 



finalement le résultat (2.1). 



Définition: Quand on a une évaluation asymptotique comme au théorème 5 (respectivement 
6) on dit que la convergence est linéaire (respectivement quadratique). 

2.5.3 Méthode de Steffensen 

Le 5 2 d'Aitken appliqué à (x n ) permet d'engendrer une suite (x' n ) qui converge plus rapidement 
que (x n )\ il est donc naturel de poursuivre les itérations en remplaçant (x n ) par (x' n ), puis calculer 
g{x' n ). On procède donc comme suit: 



Vn+1 = 9{Vn) 

Vn+2 = g(Vn+l) 



ce qui permet d'engendrer la suite: 



j _ „ (Vn+l-Vn) 

L n+1 — V-n 



Vn+2 ~ 2y n + l + Un 
x , (g«)-<) 2 

n g(g{x' n ))-2g(x' n ) + x' n 
d'où le schéma suivant connu sous le nom du schéma de Steffensen: 

xq donné 

Sn+i = G(x n ) 

„, , { % si g(g(x)) - 2g(x) + x = 
G(x) = 1 (g(x)-x) 2 

{ X g(g(x))-2g(x)+x 

Théorème7: La méthode de Steffensen est quadratique. 

Preuve: 

2d(x)(g(x) - x)(g'(x) - 1) - (g{x) - x) 2 [d'{x)\ 

(X> (g(g(x))-2g(x)+xy 
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avec d(x) = (g(g(x)) - 2g(x) + x) et d'{x) = [g'{g(x))g'(x) - 2g'(x) + 1] 
le théorème des accroissements finis implique : 

s(s(z)) = g( x ) + g'(vn)(g(x) - x) % e [x,g{x)] 

donc 

d(z) = ff(ff( a; )) - 2 g( x ) + x = g(x) + g' (rj x )(g(x) - x) - 2g(x) + x 

(g'{^)-l)(g{x)-x) 

en utilisant ce développement on a: 

%'(%) - !)(</(%) - i) - (ff'(s(^))s'(x) - 2 ff '(ï) + 1) 



G'(x) = 1 ■ 



(s'fe) - 1) 2 



si x — > s, comme ?^. G [x,p(x)] alors 77^ — > s, il est donc facile de voir que G' {s) = si </(s) 7^ 1. 
La formule de Taylor à l'ordre deux permet d'obtenir: 

e^+i = x n+1 - s = G(x n ) - G(s) = G'(s)(x n - s) + -(x n - s) 2 G"(Ç n ) 

où Ç n (z]x n , s[; comme G'(s) = on obtient alors: 



ri— +00 ei 2 



^^ = ^"( S ) 



ce qui montre que Steffenson est d'ordre deux. 

2.6 Méthode de Newton-Raphson élémentaire (1642-1727) 

Soit à résoudre f(x) = OÙ / : 1R — > H est continûment dérivable. On peut essayer de se ramener 
au cas de la méthode itérative en remarquant que si s est solution de f{%) = alors s est aussi 
solution de x = x + kf(x) quel que soit k ÉE H, d'où l'algorithme : x n +i = x n + kf(x n ). Si k est 
judicieusement choisi ça peut marcher ( il faut que |1 + kf'{x)\ < 1); dans ce cas la convergence 
sera linéaire. 

Question: Peut— on déterminer une fonction h(x) tel que: la suite engendrée par la 
fonction g(x) = x + h(x)f(x) converge quadratiquement vers s. (c— à— d f(s) = )? 

Réponse: pour avoir une convergence quadratique il faut que: g'(s) = 0. Or g'(s) = l-\-h'(s)f(s)-\- 
h(s)f'(s) du fait que f(s) = 0, pour avoir g' (s) = il suffit de prendre h(x) = — fl ) -. 

Algorithme de Newton-Raphson: 

Choisir un xq et déterminer la suite (^n). ne IN définie par: 

_ ~ f(x n ) 

f'(Xn) 

Remarque: g(x) = x — j,f{ =$- g f (s) = 0; si f( 2 > existe et est continue alors Newton-Raphson 
est quadratique. En effet, écrivons le développement limité de / à l'ordre 2: 



f(s) = = /(»„) + f(x n )( S - x n ) + l/2/« (£„)(« - x n f 
. . _ i /' 2) K-) c 2 

S - 2 f'(î n ) Ë n ' 

e„ + i 1 f"(s) 



de ce développement on tire: e n+ \ = x n+ i — s = 2 fA\ e n e ^ P ar passage à la limite on a: 



3h» e? n 2 /'(.s) 
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y=f(x) 




Figure 2.7: Principe de la méthode Newton-Raphson 



Interprétation géométrique de la méthode de N— R 

Soit xq un point de l'intervalle [a, 6] sur lequel / possède un zéro, l'équation de la tangente Tq en 
ce point est donnée par: 

~ — = / l^o) 

X — Xq 

l'intersection de Tq avec l'axe (ox) est le point x\ d'ordonnée y = 0, donc xi = xq — {}f~ \ . L'équation 
de la tangente T\ au point x\ est donnée par: 



y-f(&i) 



: /'(*i) 



/'■Mi 



l'intersection de Tj avec l'axe (ox) est le point x<i (d'ordonnée y = 0), donc x<i = x\ — {}, } \ 
On refait la même opération plusieurs fois, on engendre la suite {x n ) n çJ^ définie par: 

xq donné 



•^71+1 Xfl 






2.6.1 Convergence globale de la méthode de Newton-Raphson 

Théorème 8: soit / E C 2 [a,b] on suppose que: 

(i) /(«)/(&) < 

(ii) pour tout x ÉE [&,b] f f ($) 7^ (/ est monotone sur [a, b] ) 

(iii) pour tout x ëE [a, b] f"(x) > ou f"(x) < (/" est de signe constant sur [a, b]). 

(iv) On choisit x e [a,b] tel que f(xo)f"(xQ) > 

alors \/xq G [a, b], satisfaisant (iv), la méthode de Newton-Raphson converge vers l'unique solution 
s de f(x) = dans [a, b]. 

Preuve: Il ne peut y avoir que quatre cas de figure: 
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i5 



*/^~~~ 




b 


t ® l 

a ® J 


. 


1 © 
Vs® 


1 — " b a "- 1 



Figure 2.8: Méthode Newton-Raphson: les divers cas possible 

1 f(a) < 0, f(b) > f"(x) < sur [a, b] 

2 /(a) > 0, f(b) < f"(x) < sur [a, b] 

3 /(a) < 0, /(&) > /"(je) > sur [a, 6] 

4 /(a) > 0, f(b) < f"(x) > sur [a, b] 

en changeant / en — / on voit que 2) et 4) deviennent 3) et 1). Il suffit donc de considérer 1) et 3). 
En changeant / en — / et x en —x le cas 3) se ramène à 1). Il suffit donc de faire la démonstration 
pour le cas 1. 

1.) Dans ce cas, /(a) < 0, /(&) > j"{x) < sur [a, b]: alors f(a) < f(b) et l'on a f ^ 
sur [a,b] donc f'(x) > pour tout x G [a, &], / est donc strictement croissante. 
D'après le théorème de Bolzano (/(a) < 0, f(b) > et / continue sur [a, b] ) il existe une unique 
racine s dans [a, b] (car /' > 0) 

Pour montrer que la suite (x n ) est convergente, on montrera qu'elle est croissante majorée. En 
effet, montrons que: 

Vn > si x n < s alors x n < x n +i < s 

Si x n < s alors f(x n ) < f(s) = donc — î,Z n \ > d'où x n+ \ = x n — \jf" \ > x n (vu que 

/'(*„) > 0). 

Par application de la formule de Taylor à l'ordre 2 on a: 



f(x„+i) = f{x n ) + f'(x n )(x n+1 -x n ) + -/"(Ç„)(£„+i -5„) 2 , J„e[£„,£„ + i] 

= /(£„) + /'(£„)(-/M) + 1/"(Ç„)(£„ +1 - 5„) 2 



/"(Ufe+i -£„) 2 <0 (/"<0) 



(2.2) 



Finalement on a: f(x n +l) < 0, £„+2 = x„+l - /'(^',©) > £ n+i et donc s »»+2 > Xn+1- 
D'autre part, on appliquons la formule de Taylor au second ordre, on a: 

o = f(s) = /(*„) + (s - x n )f(x n ) + {s ~* n) /"(£,) , £„ e [â„, S] 
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Comme /" < 0, alors: 

-f(xn) < {s - x n )f(x n ) =ï - ' *\ < s- x n 

J \ x n) 
_ ~ _ /(an) - - _, _ - 

=> .T n + i < S 

Conclusion: La suite (x n ) est croissante majorée et est donc convergente. 

Choix de xq\ Dans le cas que nous considérons (/" < 0). d'après (iv) la convergence est as- 
surée pour xq tel que: xq € [a,b] et f(xo) < 0. Dans ce cas de figure on a xq < s. D'après ce qui 
précède on a: 

Xq < X\ < x-2 < ■ ■ ■ < x n < s 

la suite converge donc vers un point s tel que s = s — -lA^i =^ f(s) = 

Remarque: On se place toujours dans le cas 1 de la démonstration précédente. 

• Si /{xq) < 0: d'après ce qui précède, la suite est convergente vers s. 

• Si f{xo) > c.à.d Xq > s. xi = xq — ï>^[ < Xq, en appliquant la formule de Taylor, on a: 

/(»i) = /(£o) + ( x i ~ x )f{xQ) + -(xi - xo) 2 f"(c) , c e [$lx q ] 

= f(ï n ) _ IM f(r n ) + l ( f&o) )2 f»( c ) = l ( f@o) \2 f{c) < 

S{ ° } /W [0) + 2\f'(x ) } } {C) 2\f'(x ) ) f [C) < 

ainsi, on a trouvé un autre point d'initialisation x\ tel que f{xi)f"{xi) > 

2.7 La méthode combinée 

Soit / une fonction telle que f(a)f(b) < 0, si /' et f" gardent des signes constants sur le segment 
[a,b\. En combinant la méthode de Lagrange (x n ) et la méthode de Newton-Raphson (x n ) on 
obtient une méthode dont chaque étape permet de déterminer les valeurs par défaut et par excès 
de la racine s de l'équation f(x) = 0. Quatre cas se présentent (voir figure): 

1. f'(x) > et f"(x) < pour tout x G [a, b] 

2. f(x) < et f"(x) < pour tout x e [a, 6] 

3. f(x) > et f"(x) > pour tout x G [a, b] 

4. f(x) < et f"(x) > pour tout x G [a, b] 

Nous traitons le troisème cas seulement. Les autres cas se déduisent du premier par changement 
de / en — / ou x en —x ou les deux en même temps. 

Soit f'(x) > et f"{x) > pour tout x G [a, b], posons xq = a (initialisation de Lagrange) et 
xq = b (initialisation de Newton-Raphson), d'après ce qui précède on a: 

f(x n )(x n —x n ) 

X n + 1 — X n f(x n )-f(x n ) 

Xn+1 = X n f'( x . ) ' n = ' ' ' ' ' ' 







*K')I 
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27 



Figure 2.9: Méthode combinée. 

A chaque étape la méthode de Lagrange est appliquée au segment 1%,^]. Dans ce cas de figure, 
on a l'encadrement suivant: 

Si on souhaite avoir l'approximation de s par x n à une précision e près, on arrête les itérations dès 
que: x n — x n < e. A la fin des itérations, le mieux est de prendre comme approximation s de la 
racine la moyenne arithmétique de x n et x n c'est-à-dire: 



« — ~ y^n \ %n) 
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2.8 Execices 

Exercice 1. Résoudre dans H l'équation suivante: 1.110 x + 2x — 1.4 = 0. 

Distinguer le cas de la simple précision du cas de la double précision, conclure. 

Exercice 2. Soit / une fonction continue sur [a, b] vérifiant f(a)f(b) < 0. 

a.) Démontrer que l'algorithme de dichotomie engendre une suite {c n , n ÉE IN} telle que : 

. b — a 
\ c n — s < -r , n > 1 

où s est une solution de f(x) = 0. 

b.) Déterminer approximativement le nombre d'itérations pour résoudre f(x) = x 3 + 4x 2 — 10 = 

avec une précision de e = 10~ 5 pour a = 1 et b = 2. 

c.) Etendre le résultat de b.) pour résoudre f{x) = dans [a, b] quelconque. 

Exercice 3. Montrer que les fonctions suivantes vérifient les conditions de Lipschitz avec < L < 1: 

• f(x) = 2 + ±\x\ xe [-1,1] 
. f{x)= 1 - «€[2,3] 

• f(x) = 5 - | cos(3^) x e [0, 3f] 

Exercice 4. Soit m un nombre réel et soit e G 1R tel que: \e\ < 1. Montrer que l'équation 

x = m — e sin x 
(dite de Kepler) a une solution dans l'intervalle [m — 7r,m + tt]. Cette solution est elle unique? 

Exercice 5. Quelle est la valeur de: 

2 + J 2 + ^/2 + .7. et s = 1 

a+ -; — ^r 



Exercice 6. Tracer la courbe y = g(x) = y/Q + x et la droite d'équation y = x. déterminer les 
points fixes de g. Soit la suite récurrente définie par x n +i = g(x n ), avec le choix de xq = 7 calculer 
Xi, X2- La suite est-elle convergente? 

Exercice 7. Résoudre l'équation 5x = e x par: 
a.) la méthode des approximations successives. 
b.) la méthode de Newton. 

Exercice 8. Montrer que l'équation .t 3 + 2x — 2 n'a qu'une solution réelle, laquelle est très proche 
de 0.77. Appliquer la méthode des approximations successives à: 

a.) x n+1 = l-xl/2 ; 

2 

b.) Xj,a_-\ = — : 

1 n+1 xl + 2 ' 



Exercice 9. (LEONARD DE PlSE.) On se propose de résoudre 1' équation de Léonard de Pise (vers 
1225) donnée par: 

f(x) = x 3 + 2x 2 + 10a; - 20 = (LP) 
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a.) Montrer que l'équation (LP) admet une seule solution contenue dans [1,2]. (En réalité cette 

solution est très proche de 1.36) 

b.) Considérons les trois méthodes d'approximations successives suivantes: 

mi : x n+1 = jg(-a£ " 2.t^ + 20) , x e [1, 2] 

m 2 : x n+1 = xl + 2x\ + llar„ - 20 , x e [1, 2] 

20 . . 

■m-i : x„_i_i = — , in G 1,2 

3 7l+1 xl + 2x n + 10 l ' J 

Etudier la convergence et la limite de chacune de ces méthodes. Dans le cas où la méthode est 

convergente, préciser l'ordre. 

c.) Pourriez-vous proposer une autre méthode itérative? (On n'étudiera pas sa convergence) 

d.) Appliquer la méthode de Newton-Raphson à l'équation (LP) tout en précisant la valeur initiale 

Xo, et donner les trois premières itérations: X\, x 2 et x$. 

Exercice 10. D'après l'algorithme d'Aitken si (x n ) converge vers s tel que s — x. n ^ pour tout 
n > et s'il existe un réel \0\ < 1 tel que lim TÎ _ s . 00 ^ n+1 — Q alors la suite (x' n ) définie par: 

/ \Xn-\-\ X n ) 



Xji-\-\ ^"^71+ 1 "T" Xji 



s — x lt 

1 S — X n 



converge plus rapidement vers s dans le sens que: lim^. 

a.) Si x n = -^, montrer que x' n = pour tout n > 0. 

b.) Si x n = 1-, montrer que x' n = ,„ ] -, pour tout n > . Y-a-t-il accélération de convergence 

dans ce cas? 

c.) Si x n = (2 ,, 1 _ 1) , montrer que x' n = (4 , i+ 1 1 _ 1) . 

Exercice 11. Soit à résoudre l'équation: x — e^' = (E). a.) Montrer que l'équation (E) 

admet une seule racine s dans l'intervalle / = [|,2]. 

b.) Déterminer approximativement le nombre d'itérations qu'il faut faire pour résoudre (E) par la 

méthode de dichotomie sur l'intervalle / avec une pécision de e = 10~ 6 . 

c.) On désire appliquer la méthode des approximations successives à (E) définie par: x n +i = 

g(x n ) = e 1 -^'. 

cl.) Cette méthode est-elle convergente sur II 

c2.) Si oui, déterminer le nombre d'itérations qui permet d'assurer que 

\x n - s\ < e = HT 6 . 

d. Dans ce cas de figure, quelle est la meilleure méthode: la dichotomie ou l'itérative? 

Exercice 12. 



On se propose de résoudre f(x) = x 2 — e x = 0, x G [—1, 0]. 

1.) Montrer que / admet une solution unique s G [—1,0]. 
2.) En appliquant la méthode de dichotomie déterminer le nombre n tel que: 

\Cn~s\ <10" 4 
et calculer cq, Ci et c-j. 
3.) Soit la méthode itérative définie par: 

x n+1 =g(x n ) = -e lx/2) 
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a.) Cette méthode est elle convergente sur [—1,0]? 

b.) Si oui déterminer le nombre d'itérations à faire pour avoir une précision de 10~ 4 sur s 

(on prend xq = —0.5) 

c.) Calculer ici, x 2 et x% 

4.) Quelle est la meilleure méthode, dichotomie ou itérative? 

5.) Proposer une autre méthode itérative et étudier sa convergence. 

Exercice 13. On se propose de résoudre f(x) = Sx 2 — e x = 0, x ËE [—1, 0]. 

1.) Montrer que / admet une solution unique s ÉE [—1,0]. 

2.) En appliquant la méthode de dichotomie, déterminer le nombre n tel que: \c n — s\ < 10 et 
calculer cq, c i e "t 02- 

3.) Soit la méthode itérative définie par: 

x n+ i = g a (x n ) = ae^W 

a.) Pour quelle valeur de a la fonction g admet comme point fixe la solution s de f(x) = 0? 

On note qq cette valeur. 

b.) La méthode itérative définie par g aQ , est-elle convergente sur [—1,0]? 

c.) Si oui, déterminer le nombre d'itérations à faire pour avoir une précision de 10~ 6 sur s 

(on prend xq = —0.4) 

d.) Quelle est l'ordre de cette méthode? 

e.) Calculer x±, x- 2 et x% 

4.) Quelle est la meilleure méthode, dichotomie ou itérative? 

5.) Proposer une autre méthode itérative et étudier sa convergence. 

6.) Ecrire l'algorithme de Newton-Raphson pour résoudre f(x) = et étudier sa convergence 

sur [-1,0]. 

Exercice 14. Soit à résoudre l'équation f(x) = x — e~ x = (E). 

a.) Montrer que / admet une seule solution qu'on localisera entre deux entiers successifs p , p-\- 1 

b.) Montrer que la méthode itérative définie par: x n+ i = g(x n ) converge sur / = [p,p-\- 1]. 

c.) Montrer que la méthode de Newton-Raphson est convergente sur I. Pour votre choix de Xq, 

donner #i, x 2 et x$. Au bout de ces trois itérations quelle est la précision que vous avez atteinte 

sur l'évaluation de la racine de {E)7 

Exercice 15. Soit / : H* — * R définie par: 

/(*) = ^1 - e- 

X 

1.) Etudier la fonction / et représenter graphiquement Tf. En déduire le nombre de racines de 
f(x) = 0. Chacune des racines sera localisée entre deux entiers consécutifs. 



2.) Etudier les méthodes itératives suivantes: 



xq donne 
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a. g(x) = 1 + xe x 

b. g(x) =ïog( -) 

x — 1 

c. g(x) = (x - l)e x 

On étudiera la convergence vers des solutions du problème f(x) = en précisant le choix de 
xo et l'ordre de convergence. 

3. Ecrire la méthode de Newton-Raphson associée à l'équation f(x) = 0, choisir Xq pour que 
cette méthode converge vers la plus grande des racines. 

Exercice 16. Soit g une fonction définie sur [a, b] telle que: 

• g est dérivable sur [a, b] 

• g([a,b]) c [a, b] 

a.) Montrer que si max.\g f (x)\ x ^i at} } < 1 alors g est contractante dans [a, b]. 

b.) Dans ce qui suit, on suppose que: < L = max[</(x)| ;C( =r a y < 1, montrer que la suite (# n ) n =lM 

définie par: sc n +i = g{%n) converge vers un point s que l'on précisera. 

c.) Montrer que cette limite s est unique. 

d.) Montrer que pour tout n €= IN \x n+ i — x n \ < L n \xi — xq\. 

e.) Conclure que \x n — s\ < ïztl\ x i ~ x o\- 

Exercice 17. Soit g une fonction définie sur [a, b] telle que: 

• g est dérivable sur [a, b] et pour tout x G [a, 6], </(#) < 

• g est contractante dans [a, 6] (0 < L < 1 la constante de contraction.) 

• g([a,b]) c [a,b] 

a.) Montrer que la suite (x n ) définie par: x n+ i = g(x n ) , xq G [d,b], converge vers un point s que 

l'on précisera. 

b.) Si xq < s, montrer que pour tout n G IN on a: 

Xq < x 2 < X4 <...<... < x 2n < • • ■ < s < . . . < x 2n +i < ■ ■ ■ < X5 < X3 < Xi 



Exercice 18. On veut appliquer la méthode de Newton-Raphson (N.R) pour calculer la racine 
carrée d'un nombre réel strictement positif a. 

a.) Montrer alors que pour résoudre f(x) = x 2 — a = par la méthode de (N.R), on génère la 
suite (x n ) définie par: 

xq donnée 

-^ïï + l — 2(^)1 "T ~^~) 

b.) Montrer la convergence de la suite (x n ) si on initialise convenablement. 
c.) Montrer que la suite (x n ) définie dans b. satisfait: 

Xn-V<X ( X - ^ \2" 

— =1 — I , pour tout n G JM 

X n + V a V x + V a } 

d.) Montrer en utilisant c. que x n converge vers y/a. pour n'importe quel choix de xq. 
e.) Montrer que la convergence de x n vers s/ôt. est quadratique. 

Exercice 19. Pour évaluer y/a connaissant a(a > 0), on propose les deux méthodes suivantes: 
f xq donnée [ y donnée 

^1 < 2ax !l M? 

x\-\-a 

Pour chacune de ces méthodes 



1 T _ 2ax, l M-2 \ y nto 2 \ 
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1. Montrer la convergence vers ^/a si on initialise convenablement; 

2. Donner nn intervalle d'initialisation; 

3. Déterminer Tordre de convergence; 

4. Choisisser un xq et un yo dans le cas a = 3. 

Exercice 20. Soit / une fonction de caisses C 2 sur [a, 6], on s'intéresse à la résolution de f(x) = 

par la méthode de Newton-Raphson (N-R). On distinguera deux cas: le cas où la racine est simple 

et celui où cette racine est multiple. 

I-Racine simple: Soit s G [a, b] une racine simple de f(x) = 0. La méthode de N-R appliquée à 

f(x) = engendre la suite (x n ) définie par: x n+ i = g(x n ) avec g(x) = x — f(x)/f'(x). 

a.) Calculer g'(x), g"(x) et en déduire g'(s) et g" {s). 

b.) Quel est l'ordre de cette méthode? 

c.) On pose e n = x n — s l'erreur associée à la méthode de N-R, démontrer qu'asymptotiquement 

et au second ordre en e n on a: 

,. e„ + i _ 1 /"(s) 
-h™ el 2/'( S ) 

ll-Racine multiple: Dans cette partie, on suppose que s est une racine multiple, soit M son ordre 
de multiplicité (c.à.d que f(x) = (x — s) h(x), avec h(s) ^ ). La méthode de N-R appliquée à 
f(x) = engendre la suite (x n ) définie par: x n+ i = g(x n ) avec g(x) = x — f(x)/f'(x). 

a) Calculer g'(x) pour tout x G [a, b] et en déduire que lim^^ s </(#) = (M — 1)/M . 

b) Démontrer qu'asymptotiquement et au premier ordre en e n on a: 

hm = — — — 

n^ + oo e n Ad 

c) Montrer que s est racine simple de t(x) = f A- -, 

d) En déduire de I que la suite zc n +i = x n — J^4 converge quadratiquement vers s 

e) La méthode de Newton-Raphson modifiée est définie par: x n+ i = x n — M 1,^1 - , où M est Tordre 
de multiplicité de la racine s de /. Montrer que la méthode de Newton-raphson modifiée est d'ordre 
deux. 

Exercice 21. Soit / une fonction définie et dérivable sur [0, 1] telle que /(0) < < /(l) et 

< a < f'(x) < b ou a et b sont deux constantes. 

a.) Montrer qu'il existe s G [0, 1] unique tel que f(s) = 0. 

b.) Montrer qu'il existe une constante M telle que la méthode itérative définie par x n+ i = g(x n ) 

où g(x) = x + Mf(x) converge vers s. 

c.) Pour quelle valeur de M la suite (#n) n c]N définie ci-dessus converge quadratiqument vers s? 

Exercice 22. Soit g une fonction possédant une dérivée troisième g^> et telle que s = g{s), 
g'(s) = g" (s) = 0. Montrer que dans ce cas la limite de -^^ existe pour x$ au voisinage de s, où 

(cette convergence est appelée convergence cubique). 

Application: Soit 

. x A + bx 

cx À + a 

déterminer les relations que doivent vérifier b, c, d et s de telle sorte que x n +i = g{x n ) converge 
cubiquement vers s^0 
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Aide: 

„ _ 2x(bc - d){-3d + ex 4 ) 

9 {x > ~ (d + ex?)* 

Exercice 23. Soit / G C r ([a, b]), s une racine simple de /, on pose h(x) = [f(x)] r , r > 2. 

a. Soit g(x) = x — uihL ■> calculer g'{x) et déterminer la limite de g(x) lorsque x — > s 

b. Soit (x n ) la suite engendrée par l'algorithme de Newton-Raphson pour résoudre h(x) = 0. Quel 
est l'ordre de convergence de cette méthode?. 

c. Pour quelle valeur de m la suite X n +i = x n — m uf^K conv erge quadratiquement vers s? 

Exercice 24. La méthode de Newton-Raphson pour résoudre f{x) = perinet de générer une 
suite (x n ) définie par: 

_ f(x n ) 

Regula-Falsi propose, qu'au lieu d'évaluer f'{x n ) à chaque itération, d'utiliser une approximation 
de f'(x n ) par: 

t l, v _ f( x n) ~ f{x„-l) 
*£jl Xn— 1 

a.) Montrer que dans le cas de l'approximation de Regula-Falsi on a: 

X n -if(Xn) -X n f(x n -i) 



%n+l 



f(x n ) - f[x n -\ 



b.) Si on applique la méthode de Régula Falsi pour résoudre x 2 — 1 = 0, on supposons que 

e n = x n — 1 — » 0, montrer qu'on a: 

,. e n +i 1 

iim = — 

n^+oo e n e n -i 2 

Exercice 25. On veut appliquer la méthode de Newton-Raphson (N.R) pour calculer l'inverse 
d'un nombre réel q^O. 

a.) Montrer alors que la résolution de f(x) = ■ a = par la méthode de (N.R), on génère la 

suite (x n ) n <=]\[ définie par: 

' adonnée (23) 

b.) Montrer la convergence de la suite (#re)„ e ]N s ^ on hiitialise par un xq vérifiant: < xq < 2cx . 
c.) Si la suite (x n ) définie par 2.3 est initialisée par xq = 1, montrer que: 

_ _i-(i-«) 2 " 



d.) En déduire que la suite ainsi générer converge pour < a < 2 et que la convergence est 

quadratique. 

Exercice 26. 1.) Montrer que la fonction f(x) = 1 — 3e~ x a une racine unique s sur IR. 

2.) Par un procédé de dichotomie, trouver un intervalle de longueur 1 contenant cette racine s. 

Pour l'approcher, on propose les schémas suivants: 

3.) X e [1,2], X n + i = X n +f(x n ). 

Cette suite converge- 1- elle vers s? pourquoi? 

4.) Trouver une valeur de A G IR qui assure la convergence de la suite x n+ i = x n + Xf(x n ). 
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5.) Appliquer à f(x) la méthode de Newton-Raphson pour trouver s. 

Sur quelle fonction définit-on les itérations de point fixe? 

Dans quel intervalle doit-on choisir la valeur de départ pour que la convergence soit assurée. 

Exercice 27. Soit / une fonction de classe C A [a, b]. On suppose que / admet une racine s unique 
sur [a, 6] tel que pour tout x G [a, b]: f'(x) > et f"(x) > 0. 

1.) Donner l'interprétation géométrique de la méthode de Newton-Raphson et montrer que l'on 
génère une suite x n définie par: 

x n+1 = x n - n = g N {x n ) (2.4) 

/ { X 7lJ 

2.) Montrer qu'il existe £„ entre x n et s tel que: 

/(*„) + (s - x n )f(x n ) + { -^£f"(Xn) + ^p^/ (3) (U = (2.5) 

3.) En négligeant les termes du second et du troisième ordre dans (2.5) et en supposant que 
s « x n+ i retrouver (4.8). 

4.) En négligeant le ternie d'ordre 3 dans (2.5) et en résolvant l'équation du second degré en 
Y = (s — x n ), montrer que l'on obtient les racines suivantes: 



Y h2 = s- Xn = ^ (2.6) 

5.) Montrer qu'à partir de la solution Y\ et en supposant que s ~ #rH-ij on obtient la formule de 
Cauchy suivante: 

x n+1 = x n = gc{x n ) (2.7) 



/'(.„)[! + ^^gp] 



6.) En remplaçant »/l — " fi 7 ,) ^2 P ar 1 — f/'f )i 2 dans ^ a forixiiile de Cauchy, montrer que 
l'on retrouve la formule Halley: 

7.) Montrer qu'à partir de la formule de Halley on obtient la formule de Tchebychev suivante: 

J\ x n) r , J\ x n)j v^nj-i / \ (n c\\ 

Xn+1 = Xn ~7M {1+ 2(/'( In )) 2 ] = 9T{Xn) (2 ' 9) 

8.) Calculer g T (s) et g T (s). quel est l'ordre de la méthode de Tchebychev? 

9.) Ecrire un programme en Pascal permettant de résoudre f(x) = par la méthode de Tchebyshev. 
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